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A Bose gas in a single-beam optical dipole trap 
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We study an ultracold Bose gas in an optical dipole trap consisting of one single focused laser 
beam. An analytical expression for the corresponding density of states beyond the usual harmonic 
approximation is obtained. We are thus able to discuss the existence of a critical temperature for 
Bose-Einstein condensation and find that the phase transition must be enabled by a cutoff near 
the threshold. Moreover, we study the dynamics of evaporative cooling and observe significant 
deviations from the findings for the well-established harmonic approximation. Furthermore, we 
investigate Bose-Einstein condensates in such a trap in Thomas-Fermi approximation and determine 
analytical expressions for chemical potential, internal energy and Thomas-Fermi radii beyond the 
usual harmonic approximation. 

PACS numbers: 05.30.Jp, 37.10.De, 03.75.Hh, 51.10.+y 



I. INTRODUCTION 

Since the first observation of Bose-Einstein condensa- 
tion in dilute atomic vapors in a remarkable series of 
experiments 0, 0, 0] , ultracold quantum gases represent 
an immensely active field of research for experimental 
scientists and theorists alike. These systems offer an un- 
precedented variety of possibilities to manipulate and in- 
vestigate many-body quantum phenomena 0, [H, Q . 

A common approach to reach the Bose-Einstein con- 
densation phase transition involves laser cooling of the 
atomic vapor after which atoms in weak-field seeking 
states are transferred into a magnetic trap @, Q , followed 
by evaporative cooling. While for magnetic traps the po- 
tential depends on the magnetic substate of the atoms, 
optical dipole traps have been used to store atoms inde- 
pendently of their magnetic substate. To give a recent 
example, Bose-Einstein condensation of Strontium has 
been achieved in crossed laser beams [H, [l(| ■ 

A focused laser beam provides an almost perfect re- 
alization of a conservative trapping potential for neutral 
atoms [ll| , and is the basis for the growing field of ultra- 
cold gases in optical lattices [HI]. The first attempts to 
obtain quantum degeneracy in an optical dipole trap were 
made by Chu et. al. [H[ > wn0 employed two crossed laser 
beams. Multiple implementations of optical and evapo- 
rative cooling techniques may be used, which allow phase 
space densities close to quantum degeneracy [3, [HI, [lj| . 
Chapman et al. managed to condense 87 Rb in a dipole 
trap consisting of two crossed C02-laser beams by means 
of evaporative cooling [l7| and all-optical crossed-beam 
setups have also enabled the first observation of Bose- 
Einstein condensation of Cesium and Ytterbium atoms 
[lj|. Optical traps consisting of one single focused 
laser beam were used to produce degenerate Fermi gases 
20], Fermi gas mixtures 21], to achieve condensation of 
87 Rb [13], and they are a popular way to trap and store 
cold gases for multiple needs [H, H ^1 ES IS, IHl • 

Usually the potential of such an optical trap is approx- 
imated harmonically for low average temperatures of the 
atoms compared to the trap depth [2^, Hl| • For some ap- 



plications, however, e.g. for evaporative cooling, typical 
temperatures need not be that low. Our work is moti- 
vated by an experiment of Helm and coworkers [23| who 
use a single focused C02-laser beam to cool and store an 
ultracold cloud of 87 Rb-atoms. We aim at a more de- 
tailed understanding of the evaporation dynamics and of 
the properties of the ultracold gas in a single-beam trap. 
Crucially, our analysis is based on the true dipole trap 
potential without invoking the usual harmonic approxi- 
mation. 

Our article is organized as follows. In section |TT] we 
manage to derive the density of states of the dipole trap 
based on the usual semiclassical expression without in- 
voking further approximations. For low energies we re- 
cover the well-known expression for the harmonically ap- 
proximated potential; for energies near threshold, how- 
ever, we uncover a singularity. Subsequently, we will 
discuss the critical temperature for Bose-Einstein con- 
densation in such a trap in section IIIII Remarkably, 
without further thought the singularity in the density 
of states prevents Bose-Einstein condensation from oc- 
curring. However, a cutoff - most likely introduced by 
gravity - removes this singularity and paves the way for 
Bose-Einstein condensation. In section IIVI the influence 
of the modified density of states on evaporative cooling 
will be studied and significant differences to the usual 
harmonic approximation will be revealed. Finally, in 
section [V] we investigate the condensate wave function 
in Thomas-Fermi approximation in such single-beam op- 
tical dipole traps. We conclude with final remarks in 
section IVII 



II. DENSITY OF STATES 

We discuss an optical dipole trap consisting of a single 
focused Gaussian far detuned laser beam. T he p otential 
that affects atoms in such a trap is given by [291 ] 



Udi P (p,z) = - 



U 



(1) 
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which displays an azimuthal symmetry along the prop- 
agation (z-) axis. The energy Uq > is the maximum 
depth of the trap at p — 0, z = 0. Here, zq denotes the 
so called Rayleigh length and wq the beam waist of the 
laser w hich is c onnected to the laser wave length through 
wo = ^Azq/tt- 
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FIG. 1: Dipole trap potential (black, solid) and its harmonic 
approximation (blue, dashed), at z = y = (left) and x = 
y = (right). 

If the average thermal energy k^T of the atoms in the 
trap is much lower than Uq, the potential can be approx- 
imated harmonically which leads to the expression 



2p z 

U ha (p, z) = U [-l-\ 2 + ~2 

Wq zf 



(2) 



This harmonic approximation is usually employed to de- 
scribe cold atoms in dipole traps [2(| • 

For a proper discussion of the thermodynamic proper- 
ties of a Bose gas we aim at the true density of states g(E) 
[H of the dipole trap; the number of single-particle eigen- 
states between E and E + dE is then given by g(E)dE. 
In the usual semiclassical approximation the phase space 
is divided into Planck cells of size (2ttH) 3 such that the 
density of states becomes 
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In the following it turns out to be convenient to work 
with a scaled energy (which is zero at the trap center 
and one at the threshold) and the accompanying scaled 
density of states: 



E + U 
U ' 



g(e) = U g(E), 



(4) 



such that g(e)de — g(E)dE. 

Remarkably, inserting the full dipole trap potential ([1]) 
in the semiclassical formula @ allows us to express the 
density of states in terms of known special functions. We 
find 



gdip(e) = G 



Ve^SF ^arctan (^\J ^ " 
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Here, G 



47T 2 (2m!7o) 3/2 
9(27rri) 3 



zqWq is a dimensionless con- 



stant (see below for a more convenient expression) and 
F{<f>, m) and E((f>, m) denote the elliptic integrals [30] of 
the first and second kind, respectively, with 



F(<t>,m) 
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Expanding expression ([5]) up to second order for small 
energies leads to the well-known density of states of the 
harmonic oscillator potential ([2]), as expected. We find 



limgdip(e) 
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with the geometric mean u> of the trap frequencies de- 
rived from the harmonic approximation lp]). Thus, the 

prefactor in can also be written as G = j^jj which 
is typically a large number. The harmonic approxima- 
tion © agrees fairly well with the true density of states 
of the dipole trap potential <7dip(e) for energies e up to 
about one third of the potential depth as can be seen 
in Figure [H Clearly, for larger energies the two density 
of states differ significantly. Expression ([5|) displays a 
rj - 1 — singularity (and an additional logarithmic singular- 



ity from F(4>,m)) as e 
formula we found 



1. As a simple interpolating 
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1 (see Fig. EJ. 
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FIG. 2: Comparison of the densities of states of the dipole 
trap potential (black, solid line) and its harmonic approxima- 
tion (blue, dashed line), along with the simple approximate 
formula (dotted line). 

The singularity at the potential edge has profound im- 
plications for the thermodynamic properties of a Bose gas 



3 



in this trap. In particular, we next want to investigate its 
influence on the critical temperature Tq for Bose-Einstein 
condensation. 



tion /(e) with zero chemical potential. In other words, 
using our scaled units, the critical temperature for Bose- 
Einstein condensation for a given total atom number iV 
in the trap is here determined by [3l| . 



III. CRITICAL TEMPERATURE 

Bose-Einstein condensation occurs once the total num- 
ber of atoms exceeds the maximum possible number of 
atoms in excited states which is set by the Bose distribu- 



N= / deg(e)(cxp(eUo/k B T c )-l) 



(8) 
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FIG. 3: Comparison of the number density g(e)f(e) of the dipole trap potential (black, solid line) and its harmonic approximation 
(blue, dashed line) for r\ = Uo/k^T = 2.5,5,10,20. The areas below the graphs must be finite to allow for Bose-Einstein 
condensation. 



In case the integral diverges there is no Bose-Einstein 
condensation as for the well-known case of the one- 
dimensional harmonic oscillator or the two-dimensional 
box potential. Figure [3] shows the integrand g(e)f(e) 
of expression ((SJ) with the Bose distribution /(e) = 
(exp(rye) — 1) _1 for different ratios rj = l/o/hsT of trap 
depth to thermal energy. 

The figure reveals that the number density g(e)f(e) di- 
verges for the dipole trap near the potential edge while 
it stays finite for its harmonic approximation. Well be- 
low the singularity, i.e. for small energies, the agreement 
between the two improves for increasing 77, i.e. for lower 



temperatures: then the atoms are closer to the trap min- 
imum, where the potential is well characterized by the 
harmonic approximation. For higher temperatures the 
atoms populate states of higher energy, where the har- 
monic potential differs significantly from the dipole trap. 
For e — ► 1, the product g(e)/(e) keeps the singularity 
from the density of states and it turns out that it is not 
integrable. 

We have to conclude that - strictly speaking - the Bose 
Einstein phase transition does not occur in a single-beam 
optical dipole trap, no matter how cold the gas; which 
contradicts findings of recent experiments [22| . 
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FIG. 4: The dipole trap potential without gravity (black, 
solid), including gravity (red, dashed), and the pure gravita- 
tional potential (green, dashed-dotted) at y — z = 0. Due 
to gravity, atoms may escape from the trap as soon as their 
energy is e > 1 — Ae. Thus, the real potential is cut off at 
e t = 1 - Ae. 

In order to shed light on this curious result we have to 
look at the true experimental conditions in more detail. 
A solution of this apparent contradiction can be found 
in gravity - which we neglected up to now. Including 
gravity, the true potential reads 

TT f* U ° ( - 2 (^ + y 2 ) ^ 

Ud ^ {r) = - {i + z*/zD exp Ua(i + *V*b a ) ) ~ mgx ■ 

(9) 

As displayed in Figure [U gravity lowers the threshold en- 
ergy e t an atom needs to escape the trap in one spatial 
direction. Based on the assumption of sufficient ergodic- 
ity, we simply assume that gravity effectively leads to a 
cutoff at a certain scaled energy e t = 1 — Ae. The cut- 
off parameter Ae depends on the absolute trap depth Uq 
and on the mass of the atoms. 

As a consequence of this cutoff, the integral over the 
number density in equation ([5]) extends up to e t = 1 — Ae 
only. The integral remains finite and Bose condensation 
may take place. 

In fact, the critical temperature for a given number of 
atoms can then well be estimated as usual within the har- 
monic approximation. To give an example, in an optical 
dipole trap with Uq — 200/iK, the critical temperature 
for 10 5 87 Rb-atoms is approximately 800nK, which is in 
accordance with rj ss 250. This implies that the atoms 
are trapped so deeply that the harmonic approximation is 
clearly valid - see Figure [TJ The same conclusion can be 
drawn from Figures [3J while for fairly high temperatures 
(rj = 2.5, 5) the number density in the dipole trap differs 



from the harmonic approximation over the whole range 
of energies, this seizes to be true for low temperatures 
(77 = 20): here the only difference is the sharp singularity 
near threshold (e w 1) which, as we have argued, turns 
out to be irrelevant due to gravity. 

IV. EVAPORATIVE COOLING 

In the last section we came to the conclusion that even 
though the integral ([5]) is not integrable in the usual lim- 
its, a critical temperature for Bose-Einstein condensa- 
tion exists due to gravity or any other interaction that 
effectively cuts off or blurs the singularity as e — > 1. 
The harmonic approximation of the potential is sufficient 
to determine T c because temperatures are typically low 
enough. However, for processes involving higher temper- 
atures, as for instance evaporation dynamics, the num- 
ber density g(e)f(e) in the dipole trap differs significantly 
from the harmonic one and deviations from the harmonic 
approximation might become relevant. 

Evaporative cooling in optical traps is a simple and 
efficient method for reaching low enough temperatures 
to produce degenerate Bose [13] and Fermi [2(| gases. 
The technique is based on the preferential removal of 
atoms with higher energy than the average energy. Sub- 
sequently, the remaining atoms thermalize through elas- 
tic collisions, which leads to a lower average energy and 
temperature. Our discussion here follows [32j which de- 
scribes evaporation in a classical approximation. The 
latter is justified as typical temperatures are still much 
higher than the quantum level spacing of an atom in the 
trap potential. On the other hand the temperatures shall 
be low enough to work in the s-wave scattering regime 
with an energy independent cross section a = 8ira 2 , 
where a is the scattering length. 

A central ambition is to determine the evolution of 
the energy distribution of a trapped gas during so-called 
plain evaporation [32]. Quite generally, well above the 
transition temperature, the trapped gas is well described 
by a classical phase-space distribution /(r,p), which is 
normalized according to 

N = J d 3 rd 3 p f(r,p) . (10) 

The time evolution of /(r, p) is determined by the Boltz- 
mann equation [33j 

f £ Vr - V r C/V p + f(r, p) = J(r, p), (11) 

with the collision integral X for s-wave scattering and 
thus an energy-independent cross section a given by 



I 

X(r ' P4) = 27rm(Lfi)3 / d V^{/(r, P i)/(r, p 2 ) - /(r, p 3 )/(r, p 4 )} . (12) 
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As in [lH , P3 and P4 are the momenta of two atoms be- 
fore the collision, p! = P/2 + q' and P2 = P/2 q' the 
momenta after the collision. Here we denote a relative 
momentum with q = (p3 — p4)/2, the total momentum 
with P = P3 + P4, and |q'| = |q|. fi' denotes the direc- 
tion of q' with respect to q. 

The assumption of sufficient ergodicity leads to a phase 
space distribution as a function of the scaled single- 
particle energy e only, such that it can be written as 

/(r, p) = N JdeS ((U(r) + P 2 /2m)/U - e) /(e) , (13) 

where /(e) is the distribution of atoms as a function 



of scaled energy e which we define to be normalized to 
one. The substitution (fT3")) leads to a rigorous simplifi- 
cation of the Boltzmann equation pip . After applying 
(27rft)~ 3 / d 3 rd 3 p 5 (U(r) + p 2 /2m - e) to both sides of 
(fTTj) . the gradient terms on the left-hand side sum to 
zero, and only a term g(e)df(e)/dt survives. The right 
hand side can also be simplified, as described in [32] . As a 
general result, the following kinetic equation for the time 
evolution of the distribution function /(e) (initially nor- 
malized to one) of a gas in a trap with density of states 
g(e) is obtained: 



g(e 4 )/(e 4 ) = / deide 2 de 3 S(ei + e 2 - e 3 - e 4 ).g(min(ei, e 2 , e 3 , e 4 )){/(ei)/(e 2 ) - /(e 3 )/(e4)} 



(14) 



Here, dot denotes time derivative / = 4^ with respect to a dimensionless time r = Nq "2 ^3° t where No is the initial 
number of atoms. Recall that /(e) is normalized to one initially. Its time evolution is obtained numerically by writing 
the integral as a sum over a discretized energy scale. In Figure [5] we show evaporation dynamics for dipole trap (top) 
and for the truncated harmonic potential (bottom). Furthermore, the evolution of the number density g(e)f(e) is 
shown for the two potentials; its integral reflects the number of atoms remaining in the trap. 
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FIG. 5: Time evolution of normalized distribution function /(e) during plain evaporative cooling in a dipole trap (a) and a 
harmonic trap (c); time evolution of the number density g(e)/(e) for the dipole (b) and the harmonic trap (d). Initially, /(e) 
is chosen to correspond to a Boltzmann distribution with fceT = Uo, normalized to one. The arrows indicate increasing time; 
time steps grow exponentially. 
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We see from Figure[5]that for all times, /(e) can well be approximated by an exponential, i.e. a Boltzmann dis- 
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tribution. Apparently, during evaporation the thermal 
nature of the distribution is preserved. In the case of the 
dipole trap potential, it is only at the edge of the poten- 
tial well that /(e) differs somewhat from an exponential, 
due to the singularity in the density of states. 

Being a Boltzmann distribution, the slope of the 
graphs allows us to extract a temperature of the atoms 
remaining in the trap. 

Furthermore, as the graphs in Figure [5] a) and c) re- 
veal, the rate of change of /(e) for the dipole trap is 
smaller than for the harmonic potential. The same con- 
clusion can be drawn from parts b) and d) of the same 
figure. As the integral over g(e)f(e) equals the number of 
trapped atoms, one finds that the evaporation process in 
the truncated harmonic trap happens much faster than 
in the dipole trap. However, at the end of the evapora- 
tion less atoms are left in the dipole trap, while the final 
temperature is approximately the same. 

In order to investigate these observations in more de- 
tail, we calculate the number of atoms remaining in the 
trap by integrating over the number density g(e)/(e). 
Moreover, we determine the temperature of the remain- 
ing atoms by fitting a Boltzmann distribution to /(e), as 
a function of time. 

Before we discuss these results, however, we want to 
turn to the semianalytical solutions of equation (fT4")l . 
which are described in [321 ]. The starting point is the re- 
placement of the exact energy distribution function by a 
truncated Boltzmann distribution, which we saw is pretty 
well justified by looking at Figure 02 

/(e, t) « a{t) eM-U t/k B T(t))e{e t - e) . (15) 

The Heaviside step function ensures that f(e,t) = 
for e > e%; particles with an energy above threshold 
escape from the trap. It is worth mentioning that 
T = T(t) is not a temperature in the strict thermody- 
namic sense, since evaporation is a nonequilibrium pro- 
cess; but it surely is a convenient parameter characteriz- 
ing the essentially nonequilibrium distribution. For the 
sake of clarity we introduce the dimensionless "temper- 
ature" T = U-qT/Uq and a scaled total internal energy 
E = E/(UqNq). We keep the single-particle energy e 
scaled to Uq as well. 

Substituting the truncated Boltzmann distribution 
(p~5|) in equation ([14"]) , the integration can be performed 
easily. In order to determine the rate of change of the 
number of trapped atoms one integrates over untrapped 
energy states with e 4 = e 4 + e 2 — e 3 > e t > ei, e 2 > e 3 : 

/>oo 

N ev /N = - de 4ff (e 4 )/(e 4 ) = ~a 2 f e^V (16) 

Jet 

with 

V = r deg(e){(e t - e - f )e^ f + f e""} , (17) 
Jo 

which corresponds to the "effective volume for elastic col- 
lisions" of f32|. Note that the factor a = a(r) in equ. HU) 
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FIG. 6: Fraction of atoms N/No remaining in the trap as a 
function of scaled time t/G after starting evaporation from 
a scaled temperature T — k^T/Uo = 1. The black, solid 
line is the result for the dipole trap and the blue dashed for 
its harmonic approximation obtained by integration of the 
differential equations resulting from the truncated Boltzmann 
approximation. The squares and circles are the results from 
the fully numerical treatment for dipole and harmonic trap, 
respectively. 

has to be determined from the normalization condition 
a(T) = (N(T)/N )/tfdeg(e)e-^\ 

The evaporation of atoms leads to a loss of internal 
energy of the trapped gas. Following the elaborations in 
[32| . we find an evolution equation somewhat similar to 
(flu]) for the scaled total energy E. On the other hand, 
the rate of change of the internal energy due to the evap- 
oration of atoms is 

E = CT + n(N/N ) (18) 

where C denotes the heat capacity C — (dE/dT)^ and 
PL = (dE/d(N/N )) f = E/(N/N ) a chemical potential. 
We thus arrive at coupled differential equations which 
describe the rate of change of temperature, atom number 
and internal energy during evaporation in closed form. 
These can be solved numerically. 

Using this approximative set of differential equations, 
we want to investigate the evaporation process in the 
dipole trap based on the true density of states (O and 
compare with results for the harmonic approximation. 
Furthermore, we compare with a full numerical treat- 
ment based on the Boltzmann equation (|14[) . In Figure [S] 
we show the evolution of the number of remaining atoms 
and in Fig. their temperature, respectively, during 
evaporation. We choose an initial scaled temperature of 
T = k-gT/Uo = 1. Recall that time is dimensionless with 

t = Nq t with Nq the initial number of atoms in the 
trap (similar to 32], times are scaled to t/G). The curves 
are obtained by solving the above mentioned coupled dif- 
ferential equations based on the truncated Boltzmann ap- 
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FIG. 7: Scaled temperature T — k-sT/Uo of the remaining 
atoms as a function of scaled time t/G after starting evap- 
oration from a temperature T = 1. The black, solid line is 
the result for the dipole trap and the blue dashed line for its 
harmonic approximation obtained by integration of the dif- 
ferential equations resulting from the truncated Boltzmann 
approximation. The squares and circles result from the fully 
numerical treatment for dipole and harmonic trap, respec- 
tively. 



proximation, while the symbols follow from the numerical 
solution of equ. (fT4")) . by integrating over g(e)f(e) for the 
number of atoms and by fitting an exponential to /(e) to 
determine the temperature. 

Figures[6]and[7]confirm that the evaporation process in 
the dipole trap is indeed much slower than in a harmonic 
trap as already indicated in Figure [5] The harmonic ap- 
proximation may lead to predictions that differ by more 
than an order of magnitude. Furthermore, we see that in 
a dipole trap, more atoms get lost eventually compared to 
the harmonic trap. The final temperature, on the other 
hand, turns out to be pretty much the same. 

In the case of the dipole trap the truncated Boltzmann 
approximation does not describe the distribution func- 
tion /(e) accurately over the whole energy range. Clearly, 
we see deviations for energies near the threshold in Figure 
[5j These explain the differences between the numerical 
and the approximative results for atom number and tem- 
perature in Figures [5] and [7] 



V. A BOSE- EINSTEIN CONDENSATE IN A 
SINGLE-BEAM DIPOLE TRAP 

We here assume that a Bose-Einstein condensate has 
been created in or transferred to a single-beam optical 
dipole trap. At temperatures well below the critical tem- 
perature, the condensate wave function is well described 
by the mean field Gross-Pitaevskii equation 

(-^A + U(r)\ V>(r) + g\m\ 2 m = M V(r) (19) 



with n the chemical potential. Often, the so-called 
Thomas-Fermi approximation may be employed 0, Q 
which describes the BEC in the limit of large atom num- 
bers such that the contribution of the kinetic energy may 
be neglected with respect to the internal energy and the 
potential energy. The condensate wave function is then 
given by the simple formula 



|V(r)| s 



U(r) 



(20) 



valid for all positions r such that the right hand side 
is positive, and tp(r) =0 otherwise. The parameter g 
describes the strength of the atom-atom interaction and 
is connected to the s-wave scattering length a through 
g = 4irh 2 a/m. For a harmonic approximation, the 
Thomas-Fermi wave function (|20|) leads to the analyti- 
cal expression [i = (15g./V/47ru;oZo^o) 2 ^ 5 for the chemical 
potential, where \x is normalized to Uq and positive (sim- 
ilar to e) . Remarkably, in Thomas- Fermi approximation, 
even the full dipole potential |T]) allows for an analyti- 
cal treatment. By integrating expression (|2"0|) over r and 
equating it with the total atom number N we find the 
simple expression 



Air 
9fJ 



z w Uq I 



1-/X 



(4/x - 3) 



— 3(/i — 1) arctan 




(21) 



Similar to the density of states ([5]), the number of parti- 
cles N as a function of the chemical potential /i develops 
a singularity at \x = 1 in the dipole trap. Clearly, in 
the harmonic trap the particle number for any \x remains 
finite. For small condensate numbers N the chemical 
potentials agree, as shown in Figure [5] 
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FIG. 8: Chemical potential versus atom number of a conden- 
sate in a dipole trap (black, solid) and a harmonic trap (blue, 
dashed) in Thomas-Fermi approximation. 

Furthermore, we found an analytical expression for the 
internal (interaction) energy in Thomas-Fermi approx- 
imation for a condensate. By integrating the squared 
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expression (|20[) , one finds for the harmonic approxima- 

,2/5 

tion the well-known Ei nt (N) 



15g 



N 7/5 



4irz w*Uo / 

while for the dipole trap the result can only be given 
analytically as a function of the chemical potential: 



4 5 



/i 



¥ (m-i) + ¥ (m-i) 2 



+ (1 (fj, — l) 2 ) arctan 

3 



1 



(22) 



Clearly, the internal energy remains finite for any finite 
number of atoms. 

Finally, we determine the Thomas-Fermi radii of the 
condensate in such an optical dipole trap. They are de- 
fined by means of the condition U(r) = fx and character- 
ize the extensions of the condensate in the various direc- 
tions. For the harmonic potential they are well known 
to be R p {n) = "^fv^ an d ^z(aO = Z <i\/T L - For the dipole 
trap, we find easily 



ln(l-/i) ; R z (fi) = z Q 



1-M 



.(23) 



Clearly, as /i — > 0, these expressions coincide with the 
harmonic case. In general, however, in particular along 
the z-axis (direction of propagation of the laser), the 
Thomas-Fermi radius in the dipole trap differs strongly 
from the one obtained in the harmonic trap for large /i. 
Both expressions diverge for fi — > 1. In Figure [H] we show 
R z (fi) for both traps. 




FIG. 9: Thomas-Fermi radius along the propagation direction 
of the laser versus chemical potential for dipole trap (black, 
solid) and its harmonic approximation (blue, dashed). 

It should be emphasized that temperatures and atom 
numbers that may be achieved with evaporative cool- 



ing as described in the first part of this paper lead to 
"small" condensates that may well be described with the 
harmonic approximation. Nevertheless, it is conceivable 
that much larger condensates are created in some other 
trap before being loaded into a single-beam optical dipole 
trap such that deviations from the harmonic trap become 
relevant and our results apply. 

VI. CONCLUSION 

We found an analytical expression for the density of 
states of an optical dipole trap potential, consisting of 
one single focused laser beam. The density of states has 
a singularity at the potential edge of the trap. 

Without further consideration, this singularity would 
prevent Bose-Einstein condensation from occurring. 
However, since an effective cutoff will be present in 
current-day experiments (e.g. due to gravity), condensa- 
tion can be achieved. Then, near condensation, the usual 
harmonic approximation is valid and may be used to de- 
termine, e.g., the critical temperature for Bose-Einstein 
condensation. 

Still, the existence of this singularity makes itself felt 
whenever states of higher energy are involved as, for 
instance, in evaporative cooling. We found that plain 
evaporative cooling is much slower (about one order of 
magnitude) in the dipole trap compared to its harmonic 
approximation. In addition, for the optical dipole trap 
more atoms get lost on the way to the same temperature 
compared to the harmonic trap. 

Finally, we determine analytical expressions for chem- 
ical potential, internal energy and Thomas-Fermi radii 
for the true optical dipole trap in Thomas-Fermi approx- 
imation. It turns out that the chemical potential and 
the radii of the condensate develop singularities similar 
to the density of states. 

While we here concentrate on bosonic gases, it will 
be interesting to see how the semiclassical expression ([5]) 
will be of relevance for the behavior of fcrmionic gases in 
single-beam dipole traps. 
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